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Abstract 

We present some aspects of the so-called additive coalescence, with a focus 
on its connections with random trees, Brownian excursion, certain bridges 
with exchangeable increments. Levy processes, and sticky particle systems. 
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1. Additive coalescence in the finite setting 

The additive coalescence is a simple Markovian model for random aggregation 
that arises for instance in the study of droplet formation in clouds [131 [TT] , gravita- 
tional clustering in the universe jJU], phase transition for parking ^Uj, ... As long 
as only finitely many clusters are involved, it can be described as follows. A typical 
configuration is a finite sequence a:i > . . . > x„ > with Xi — 1, which may be 
thought of as the ranked sequence of masses of clusters in a universe with unit total 
mass. Each pair of clusters, say with masses x and y, merges as a single cluster with 
mass X -|- y at rate K{x, y) — x + y, independently of the other pairs in the system. 
This means that to each pair of indices with 1 < i < j < n, we associate 

an exponential variable e{i,j) with parameter Xi + Xj, such that to different pairs 
correspond independent variables. If the minimum 71 := mini<j<;j<„ e{i,j) of these 
variables is reached, say for the pair (ia,jo), i.e. 71 = e(io,jo)j then at time 71, we 
replace the clusters with labels iq and jo by a single cluster with 
Then the system keeps evolving with the same dynamics until it is reduced to a 
single cluster. 

An additive coalescent {X{t),t > 0) started from a finite number n of masses 
is a Markov chain in continuous times for which the sequence of jump times 71 < 
. . . < 7,1-1 has a simple structure. Specifically, the increments between consecutive 
coalescence times, 71, 72 — 71, ... , jn-i — ln-2, are independent exponential vari- 
ables with parameters n — l,n — 2,...,1, and are independent of the state chain 

*Laboratoire de Probabilites et Modeles Aleatoires and Institut univcrsitairc de France, Uni- 
versite Paris 6, 175, rue du Chevaleret, F-75013 Paris, France. E-mail: jbe@ccr.jussieu.fr 



16 



Jean Bcrtoin 



(X(7fe), k — 0, . . . ,n ~ 1). This elementary observation enables one to focus on the 
chain of states, and to derive simple relations with other random processes, as we 
shall now see. 

Pitman TT^ pointed at the following connection with random trees. Pick a tree 
r*^"-' at random, uniformly amongst the n"~^ trees on n labelled vertices. Enumerate 
its n — 1 edges at random and decide that they are all closed at the initial time. 
At time k = l,...,n — 1, open the fc-th edge, and say that two vertices belong 
to the same sub-tree if all the edges of the path connecting those two vertices are 
open at time k. If we denote by the ranked sizes of sub-trees at time k, 

then the chain y^"^ = fc = 0, . . .n — 1) has the same law as the chain of 

states (X("^(7fc), A; = 0, . . . ?i — 1) of the additive coalescent started from the initial 
configuration (1/n, . . . , 1/n). 




Forest derived at time 10 from a tree with 14 vertices 



Second, we lift from |S] a different construction which is closely related to hash- 
ing with linear probing, cf. jl(J|. We view the initial configuration xi > . . . > a;„ > 
as the ranked jumps of some bridge b — {b{u), < u < 1) with exchangeable incre- 
ments. That is we introduce J7i, . . . , n independent and uniformly distributed 
variables and define 

n 

- ^a:,(l{„>C/,}-u) , 0<u<l. (1.1) 




bridge b 



Next, we consider a path transformation (see the picture below) that has been 
introduced by Takacs [21] and used by Vervaat j^Sj to change a Brownian bridge 
into a normalized Brownian excursion. Specifically, we set 

e{u) = b{u + fi[modl])-b{fi-), 0<u<l, (1.2) 
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excursion e 



Finally, for every t > 0, call t-interval any maximal interval [a, 6[C [0, 1] on 
which 

tu — e{u) < max {(to — e{v))^ ,0 < v < u} , for all u G [a, b] . 




graph of V ^ tv — e(w) and t-intervals (hatched) 



It is easy to see that the t-intervals get finer as t increases and tend to reduce 
to the jump times of e when t ^ oo. Denote by F{t) the ranked sequence of the 
sums of the jumps made by e on each t-interval, and by < (5i < . . . < Sn-i the 
jump times of F{-). Then the chain {F{Sn-k-i), fc = 0, . . . n — 1) has the same law 
as the chain of states {X{jk), k — 0, . . .n — 1) of the additive coalescent started 
from the initial configuration {xi,. . . , a;„). 

2. Standard and other eternal coalescents 

Dealing with a finite number of clusters may be useful to give a simple descrip- 
tion of the dynamics, however it is a rather inconvenient restriction in practice. In 



18 



Jean Bcrtoin 



fact, it is much more natural to work with the infinite simplex 

f oo 

= <. X = {xi,X2, ■ ■ ■) : Xi > and 2, 2:^ = 1 



endowed with the uniform distance. In this direction, Evans and Pitman |12j have 
shown that the semigroup of the additive coalescence enjoys the Feller property on 

Approximating a general configuration a; S 5-'- by configurations with a finite 
number of clusters then enables us to view the additive coalescence as a Markovian 
evolution on It is interesting in this setting to consider asymptotics when the 
coalescent starts with a large number of small clusters, which we shall now discuss. 

Evans and Pitman |12j have first observed that the so-called standard additive 
coalescent {X^°°\t), —oo < t < oo) arises at the limit as n ^ oo of the additive 
coalescent process {X^^\t),—^logn < t < oo) started at time — ^logn with n 
clusters, each with mass 1/n. This limit theorem is perhaps better understood if 
we recall the connection with the uniform random tree r^"^ on n vertices that was 
presented in the previous section. Indeed, if one puts a mass 1/n at each vertex 
and let each edge have length n~^/^, then r^"^ converges weakly as n ^ oo to the 
so-called continuum random tree t^°°^; see Aldous jj. More precisely, t^°°^ is a 
compact metric space endowed with a probability measure (arising as the limit of 
the masses on vertices) which is concentrated on the leaves of the tree, and a skeleton 
equipped with a length measure which is used to define the distance between leaves. 
This suggests that the standard additive coalescent might be constructed as follows: 
as time passes, one creates a continuum random forest by logging the continuum 
random tree along its skeleton and consider the ranked sequence of masses of the 
subtrees. This yields a fragmentation process, and the standard additive coalescent 
is finally obtained by time-reversing this fragmentation process. Aldous and Pitman 
have made this construction rigorous; more precisely they showed that the tree 
T^°°^ has to be cut at points that appear according to a Poisson point process on 
the skeleton with intensity given by the length measure. This representation yields 
a number of explicit statistics for the standard additive coalescent. For instance, for 
every i g R, the distribution of x'f"^ is given by that of the ranked sequence > 
^2 > • • • of the atoms of a Poisson measure on ]0, oo[ with intensity e~*(27rx'^)~^/^(ix 
and conditioned by -f- • • • = 1. 

The continuum random tree bears remarkable connections with the Brownian 
excursion (cf. for instance Le Gall |14|'). and one naturally expects that the stan- 
dard additive coalescent could also be constructed from the latter. This is indeed 
feasible (see and also ^HI) although its docs not seem obvious to relate the fol- 
lowing construction with that based on the continuum random tree. Specifically, let 
(e(s), < s < 1) be a Brownian excursion with miit duration, and for every t > 0, 
consider the random open set 



Then G{t) decreases as t increases, and if we denote by F{t) the ranked sequence of 
the lengths of its intervals components (which of course are related to the so-called 




(2.1) 
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t-intervals of the preceding section), then (F(e~*),— oo < t < oo) is a standard 
additive coalescent. 

More generally, Aldous and Pitman 4, have characterized all the processes that 
may arise as the limit of additive coalescents started with a large number of small 
clusters. They are referred to as eternal additive coalescents as these processes are 
indexed by times in ] — oo, oo[. They can be constructed by the same procedure as 
in |5] after replacing the continuum random tree r(°°^ by a so-called inhomogeneous 
continuum random tree. 

An alternative construction was proposed in ^ and 15_. Specifically, one may 
replace the standard Brownian excursion e by that obtained by the Takacs-Vervaat 
transformation p.2|l where (&(s),0 < s < 1) is now a bridge with exchangeable 
increments, no positive jumps and infinite variation (which arises as the limit of 
elementary bridges of the type Hl.lll . see Kallenberg ^Hl)- The ranked sequence 
F{t) of the lengths of the interval components of G{t) defined by 1)2.111 then yields 
a fragmentation process, and by time-reversal, (F(e~*), — oo < t < oo) is an eternal 
additive coalescent. 

Roughly, this construction can be viewed as the limit of that presented in 
Section 1 when the additive coalescent starts from a finite number of clusters. 

3. Eternal coagulation and certain Levy processes 

A long time before the notion of stochastic coalescence was introduced, Smolu- 
chowski j2()l considered a family of differential equations to model the evolution in 
the hydrodynamic limit of a particle system in which particles coagulate pairwise 
as time passes. It bears natural connections with the stochastic coalescence; we 
refer to the survey by Aldous for detailed explanations, physical motivations, 
references ... Typically, we are given a symmetric kernel K :]0, oo[x]0, oo[— > [0, oo[ 
that specifies the rate at which two particles coagulate as a function of their masses. 
Here, we take of course K(x,y) = x + y. If we represent the density of particles 
with mass dx at time i by a measure /it(dx) on ]0, oo[, then 



where / a test function and {fJ.t,f) — J f{x)^t{<ix). Motivated by the preceding 
section, we are interested in the eternal solutions of 1)3. 1(1 . in the sense that the time 
parameter t is real (possibly negative) . It is proven in |^ that every eternal solution 
(/it)tGR subject to the normalizing condition J xfj,t{dx) = 1 (i.e. the total mass of 
the system is 1), can be constructed as follows. 
First, define the function 




{f{^ + v) - f{^) - f{y)) {x + v)^lt{<ix)^lt{dy) , (3.1) 




(3.2) 



where cr^ > and A is a measure on ]0, oo[ with J{x /\ a::^)A(da:;) < oo. We further 
impose that either cr^ > or J xA{dx) — oo. Next, let <&(•, s) be the inverse of the 
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bijection q „2 j^[sq) + q. One can check that $((7,6*) can be expressed in the 
form 

Hq,e') ^ f (l-e-9-Vt(dx), g>0, (3.3) 

J]0,oo[ 

where {pt)t£'R is then an eternal solution to Smoluchowski's coagulation equation. 
For instance, when ^ 1 and A = 0, ^ is a standard Brownian motion and we 
recover the well-known solution 

e~* / x(i~'^^ \ 

fj,t{d.x) ~ —r= exp — ] dx , i e R, a; > . 

\'2'Kx^ V 2 / 

This invites a probabilistic interpretation. Indeed, (|3.2|l is a special kind of 
Levy-Khintchine formula; see section VII. 1 in ISI. More precisely, there exists a 
Levy process with no positive jumps, — (Cr, > 0), such that 

E (exp (q^r)) = exp (r*^2^A(g)) , q>0. 

It is then well-known (e.g. Theorem VII. 1 in 5 ) that the first passage process 

rj") inf {r>0:s^r+r>x}, x>0 

is a subordinator with 

E(exp(-<7ri^))) = cxp{-x<^{q,s)) , q,x>0, 

where the Laplace exponent <&(-,s) is the inverse bijection of q ^a^.Aisq) + q- 
Thus (|3.3|l can be interpreted as the Levy-Khintchine formula for $(-,s), and we 
conclude that the eternal solution fit can be identified as the Levy measure of the 
subordinator T*^'*) for s = e*. 

This probabilistic interpretation also points at a simple random model for 

aggregation of intervals. Indeed, the closed rang 

induces a partition of [0, 00 [ into a family of random disjoint open intervals, namely 
the interval components of G(s) — [0,oo[\T'-'^\ We now make the key observation 
that 

q-is) c t("') for < s' < s, (3.4) 

because an instant at which r s£,r + f reaches a new maximum is always also 
an instant at which r s'^r + r reaches a new maximum. Roughly, H3.4|l means 
that the random partitions get coarser as the parameter s increases; and therefore 
they induce a process in which intervals aggregate. The latter is closely related to 
a special class of eternal additive coalescents, and has been studied in [7|E] in the 
Brownian case, and in |15| in the general case. 

4. Sticky particle systems 

Sticky particle systems evolve according to the dynamics of completely inelastic 
collisions with conservation of mass and momentum, which are also known as the 
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dynamics of ballistic aggregation. This means that the velocity of particles only 
changes in case of collision, and in that heavier cluster merges at the location 

of the shock with mass and momentum given by the sum of the masses and momenta 
of the clusters involved. This has been proposed as a model for the formation of large 
scale structures in the universe; see the survey article [221 ■ We now have two quite 
different dynamics for clustering: on the one hand the ballistic aggregation which 
is deterministic, and on the other hand the additive coalescence which is random 
and may appear much more elementary as it does not take into account significant 
physical parameters such as distances between clusters and the relative velocities. 
Nonetheless, there is a striking connection between the two when randomness is 
introduced in the deterministic model, as we shall now see. 

We henceforth focus on dimension one and assume that at the initial time, 
particles are infinitesimal (i.e. fluid) and uniformly distributed on the line. The 
evolution of the sticky particle system can then be completely analyzed in terms of 
the entropy solution to a single PDE, the transport equation 



Here it(x, t) represents the velocity of the particle located at x at time t, and the 
entropy condition imposes that for every fixed t > 0, the function u{-,t) has only 
discontinuities of the first kind and no positive jumps (the latter restriction accounts 
for the total inelasticity of collisions). Provided that the initial velocity u{-,0) 
satisfies some very mild hypothesis on its rate of growth, there is a unique weak 
solution to the equation H4.1|l which fulfills the entropy condition, and which can 
be given explicitly in terms of u(-, 0). 

We assume that the initial velocities in the particle system are random, and 
more precisely 



where ^ denotes the Levy process with no positive jumps which was used in the 
preceding section. 

Roughly, the dynamics of sticky particles are not only deterministic, but also 
induce a loss of information as time goes by, in the sense that the initial state of 
the system entirely determines the state at time t > 0, but cannot be completely 
recovered from the latter. In this direction, let us observe the system at some fixed 
time t > 0, i.e. we know the locations, masses and velocities of the clusters at this 
time. Let us pick a cluster located in [0, oo[, using for this only the information 
available at time t (for instance, we may choose the heaviest cluster located at 
time t in [0,1]). We shall work conditionally on the mass of this cluster, and for 
simplicity, let us assume it has unit mass. For every r e]0, i[, denote by M{r) = 
(mi(r), m2(r), . . .) the ranked sequence of masses of clusters at time r which, by 
time t have aggregated to form the cluster we picked, so M(r) can be viewed as a 
random variable with values in 5^ . Then the time-changed processes 



dtu + udxU = . 



(4.1) 



u{r, 0) = for r < and {u{r, 0), r 



>0) ^ (er,r>0). 




— oo < s < oo 
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is an eternal additive coalescent. This was established in |H] in the case of Brownian 
initial velocity; and the recent developments on eternal additive coalescents made 
in ^ |S1 El show that the arguments also applies for Levy type initial velocities. 
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